Snake states in graphene quantum dots in the presence of a p-n junction 
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We investigate the magnetic interface states of graphene quantum dots that contain p-n junc- 
tions. Within a tight-binding approach, we consider rectangular quantum dots in the presence of a 
perpendicular magnetic field containing p-n, as well as p-n-p and n-p-n junctions. The results show 
the interplay between the edge states associated with the zigzag terminations of the sample and the 
snake states that arise at the p-n junction, due to the overlap between electron and hole states at 
the potential interface. Remarkable localized states are found at the crossing of the p-n junction 
with the zigzag edge having a dumb-bell shaped electron distribution. The results are presented as 
function of the junction parameters and the applied magnetic flux. 

PACS numbers: 73.21.La. 73.22. Pr, 73.40.-c 



I. INTRODUCTION 

The study of graphene, a single layer of hexagonal car- 
bon, has led to the discovery of new phenomena that 
highlight the unusual electronic properties of this 2D 
systenPl. In particular, the linear gapless electronic spec- 
trum, together with the chirality of carriers in this sys- 
tem is predicted to allow perfect transmission through 
potential barriers (Klein paradox^. This transmission 
has a directional character and is caused by the over- 
lap between electron and hole states across the potential 
barricF^. The effect has been investigated experimen- 
tally in p-n junctions of gated graphene samplea^HH. 

In the presence of an external magnetic field, the 
electron-hole overlap at the potential barrier (or p-n 
junction) causes the appearanc e of states that propa- 
gate along the junction intcrfacd^l^. These are known 
as snake states since, in a semiclassical view, they arise 
through the coupling between counter-circling cyclotron 
orbits on either side of the p-n junction. They may also 
arise due to the presence of inhomogeneous field^and m 
warped and folded graphen^i^EU The presence of snake 
states was found to influence the electronic properties of 
graphene-based samples in the quantum Hall regim^^^. 
Moreover, the coupling of snake states have been pre- 
dicted to modify electrical current transport near the in- 
terfaces of narrow p-n-p junction^S!_ Recently, experi- 
mental evidence was provided of the chaotic coupling of 
snake states in quantum point contactJ^^. 

In this paper we investigate theoretically the interplay 
between edge and snake states of p-n, p-n-p, and n-p- 
n junctions imposed on graphene quantum dots (GQDs). 
We study the character of the different confined states by 
looking at the probability densities. The electron proba- 
bility density can be linked to the local density of states 
(LDOS) which is a quantity that can be measured exper- 
imentally using scanning tunneling microscopy (STM). 
Measurement of the LDOS allows the probing of the 
spatial structure of the confined energy levels. Such 
measurements were recently reported for graphene quan- 
tum dot^i^. We consider GQDs created by cutting a 
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FIG. 1: (Color online) (a) Number of nearly zero energy states 
within the interval \E\ < 0.1 meV as function of armchair 
(Na) and zigzag (A^^) edge atoms, (b) Number of zero energies 
as function of Na for iV^ = 5, 10, 20. 



larger graphene sample in order to obtain electronic 
confinement in a nanometer-scale structure with well- 
defined edges. The properties of the confined states 
of such G QDs in a magnetic field have been studied 
theoreticalljEE^ as well as experimentalljff^. Note that 
such p-n junctions (but of irregular shape) are also natu- 
rally present in graphene samples when the Fermi energy 
is located around the Dirac point. They are generally 
known as puddles and have investigated with scanning 
tunneling microscopy (STM)^*^ In our calculations 
we neglect disorder which for the considered small sized 
dots will be of secondary importance. 
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(b) E=0.273 eV, Ub=0 
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FIG. 2: (Color online) Energy levels of a system with a p-n 
junction parallel to the zigzag edges in a rectangular GQD 
with Lx = 7.52 nm, Ly = 8.6 nm and ^c/^o = as func- 
tion of applied voltage Ub- The system is illustrated in the 
inset. The different colors represent the different gate volt- 
ages i+Ui, for p-type and —Ub for n-type). The red dashed 
lines correspond to the energy states given by Eq. (3) and the 
green dashed lines represent the energy levels obtained from 
Eq. (4). 



One important aspect of such graphene-based struc- 
tures is the possible existence of edge states, for which 
the wavefunctions are locahzed at zigzag terminations 
of the sample-^™^™l These states have been recently ob- 
served by STA42S129|_ 

The presence of edge states can be 
especiaUy relevant for nanometer-scale graphene struc- 
tures. In particular, depending on the geometry of the 
GQDs, the edge states can correspond to the ground state 
of the systenpZl. For GQDs of general shape, Wimmer 
et al. have shown that the edge states tend to form a 
narrow band and are generally robust with regards to 
perturbation^^. In the present case we consider the ef- 
fect of a position-dependent potential profile and an ex- 
ternal perpendicular magnetic field on the energy spec- 
trum of rectangular GQDs in the context of the nearest- 
neighbor tight-binding model. The presence of the poten- 
tial interface thus introduces additional localized states, 
i.e. snake states, which can hybridize with the conven- 
tional zigzag edge states. 

The paper is organized as follows: Section II gives a 
description of the model. In Sec. Ill we show and discuss 
the analytical and numerical results. Our summary and 
conclusions are presented in Sec. IV. 
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FIG. 3: (Color online) Probability densities corresponding to 
the points indicated by (a-h) in Fig. [2] The black vertical 
line indicates the position of the p-n junction. 



II. MODEL 

The nearest-neighbor tight-binding Hamiltonian of the 
TT electrons in the honeycomb graphene lattice can be 
written as 



n 



^{tlmclcm + h.c), 



(1) 



where is the on-site energy, is the nearest-neighbor 
coupling parameter and Cj (cj) is the annihilation (cre- 
ation) operator of the electron at a site with label j. The 
external magnetic field introduces the Peierls phase in the 
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FIG. 4: (Color online) Energy levels of a rectangular GQD 
with a p-n junction parallel to the armchair edges (see the in- 
set of the figure) as function of applied electrostatic potential 
Ub for the same parameters as Fig. [2] 
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FIG. 5: (Color online) Probability densities corresponding to 
the points indicated by (a-f) in Fig. |4] In panel (a) iJ(g) = 
indicates the eighth degenerate zero energy. The black hori- 
zontal line indicates the position of the p-n junction. 



coupling term tim = iexp (27rz$;_,„), where t is the zero- 
magnetic field coupling parameter, ^i^m = (l/^o) /J^" A- 
dr, and $0 = h/e is the magnetic quantum flux and A 
is the vector potential. For graphene one has t = 2.7 eV. 
The field is given by B = Bz and we choose the Lan- 
dau gauge as A = {0,Bx,0). Then, the Peierls phase 
for a transition between two sites I and m is = 

in the x direction and ^i.m = ±(2^/3a)$c/*I'o along the 



±y direction, where $c = •iV5a^B/2 is the magnetic flux 
threading one carbon hexagon with a = 0.142 nm being 
the C-C distance. The p-n, p-n-p or n-p-n junctions are 
modeled by assuming a position-dependent on-site en- 
ergy em — e.{m). Throughout this paper we assign the 
values e = Ub {e = —Ub) for the p (n) regions, whereas 
at the interfaces between these regions the potential is 
assumed to vary abruptly. This assumption is expected 
not to influence the results qualitatively. 



III. P-N JUNCTION: ZERO MAGNETIC FIELD 

We consider an almost square quantum dot because 
it allows us to investigate the effect of both armchair 
and zigzag edges in the same sample. We are interested 
to learn how the conflned states are influenced by the 
relative orientation of the p-n interface with respect to 
the specific type of edges. Here, the length of the rect- 
angular GQD which is terminated by armchair edges is 
defined as — [{3Na/2) — l]a and the length termi- 
nated at the zigzag edges is Ly = N^VSa where Na and 
Nz are the number of C-atoms, respectively at the arm- 
chair and zigzag edges. The total number of C-atoms in 
the rectangular GQD is iV = Na{2N^ + 1). We should 
notice that the energy spectrum of a rectangular GQD 
exhibits zero energy states^^ which are confined at the 
zigzag edges. The number of zero energy states in a rect- 
angular dot depends to the number of both armchair and 
zigzag atomj2i!. Figure [ija) shows the number of states 
with nearly zero energies (we took the number of states 
within the energy interval \E\ < 0.1 meV) as function of 
Nz and N^. Our results show that for a fixed number 
of zigzag edge atoms the number of zero-energy states 
increases with increasing the armchair edge atoms (see 
Fig. [ijb)). Notice that the number of zero-energy states 
can not exceed 2Nz- 

Now we solve the Hamiltonian (1) for a system in which 
a p-n junction is parallel to the zigzag edges of a rect- 
angular GQD (sec the inset of Fig. [2] where different 
colors represent the p-type and n-type regions which are 
respectively subjected to +Ub and —Ub gate voltages). 
For numerical purposes we take as an example Na = 36 
{Lx = 7.52 nm) and Nz = 35 {Ly = 8.6 nm) in all 
the results of this paper. The energy levels of this sys- 
tem are shown as function of the gate voltage Ub in Fig. 
[2] for zero magnetic flux ^c/^o = 0. In the presence 
of a p-n junction parallel to the zigzag edges the zero 
energy-degenerate states split into two groups of degen- 
erate states with energy: i) E = +Ub and ii) E = —Ub 
where the number of the states in each group is equal. 
Note that the energy spectrum in Fig. [2] exhibits a group 
of states which on average are almost independent of Ub- 
Figure [3] shows probability density plots for the states 
indicated by the letters (a-h) in Fig. [2] The probabil- 
ity densities in Fig. |3][a) and Fig. [sj^c) exhibit a nodal 
character across the zigzag edges and consequently the 
Ub energy shift from the p- and n-regions cancel out. 
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(a) E=0.275 eV, Ujj=0 



(b) E=0.03 eV, U. =0.4 eV 






FIG. 6: (Color online) Upper panel: The same as Fig. [2] but 
with non-zero magnetic flux $c/4'o ~ 0.1. The red dashed 
lines are the LLs of an inflnite graphene sheet which are 
shifted up (down) by Ub{—Ub)- Lower panels: Probability 
densities corresponding to the points indicated by (a-d) in the 
upper panel. The black vertical line indicates the position of 
the p-n junction. 



These levels are similar to confined states in a zigzag 
nanoribbon. The energy levels of a zigzag nanoribbon 
are described, using the continuum model, by the tran- 
scendental equatiorP^ 



ky -j— K/ 



(2) 



{e/hvpY with vp — 10^ m/s being 
the Fermi velocity. In the low energy limit we take ky = 
where Eq. (2) becomes exp{±2iLxe/hvF) = —1 and 
results in 



± 



irhvp 



0,1,2, 



(3) 



The three first electronic levels of the above relation are 
shown in Fig. [2] by red dashed lines which coincide rea- 
sonable well with the position of the constant energy lev- 
els in the spectrum. Note that the agreement is better 
for low energy where the continuum model is more ac- 
curate. For the levels where the wavefunction is spread 
out inside the dot the energy levels are approximately 
linear with Ut- As seen in Figs. [3]jb,d) the probabil- 
ity density corresponding to these levels shows an os- 
cillatory behavior along the y-direction which is due to 
the confinement by the armchair edges. For armchair 
nanoribbons the wave vector ky satisfies the condition 
ky = {uyTT / Ly) + {2Tr /3\/3a) where Uy being an integei'^. 
Using Eq. (3) we take k^ = TT{nx + 1/2) /L^ along the 
x-direction. Thus the corresponding energies in the pres- 
ence of ±Ub are proportional to ±{Ub ± hvp^k"^ + 
which results in 



± 



Ub^hvpt 



2tt 
3\/3a 



^2 ^ ^ 7r(n, 



L. 



1/2) ^2 



(4) 

These electronic levels described by Eq. (4) which are 
shown by the green dashed lines in Fig. [2] for —26 < 
Uy < —20 and — I. The above arguments describe 
reasonably well qualitatively most of the energy levels 
that are found in the numerical spectrum depicted in 
Fig. [2j Because of the finite boundaries those levels 
may interact leading to anti-crossings. Aside from anti- 
crossings, the lines describe rather well the low energy 
levels in the spectrum that decrease linearly with \Ub\- 
Figs. |3][e,f) show the electronic density corresponding 
to the lowest paired levels for Ub = 0.38 eV where the 
electrons are only confined in the n region. Figures|3jg,h) 
show those states that are influenced by both zigzag and 
armchair edges. 

The energy levels of a rectangular GQD subjected to 
a p-n junction parallel to the armchair edges is shown in 
Fig. |4j The system is depicted in the inset of Fig. |4] 
Since the p-n interface is now located perpendicular to 
the direction of the edge states (i.e. zigzag edges) the 
energy spectrum exhibits a complex behavior as func- 
tion of Ub- In this case the spectrum is not symmetric 
under switching Ub — > —Ub for \E\ < Ub which is due 
to the fact that the number of p-type and n-type atoms 
are unequal. The probability density corresponding to 
the points indicated by arrows are shown in Fig. [5] For 
Ub = and E — the carriers are confined at the zigzag 
edges (see Fig. [sja) and the level is eighteen fold degener- 
ate) . Notice that the rectangular GQD with Na — 36 and 
Nz = 35 has 16 zero energy states (see Fig. [ija)). The 
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FIG. 7: (Color online) Current density pro- 
file corresponding to the states indicated by 
(a-d) in Fig. [6] The black line indicates the 
position of the p-n junction. The inset in 
(b) shows a streamline plot of the enlarge- 
ment region. 
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FIG. 8: (Color online) The same as Fig. [Ijbut in the presence 
of magnetic flux $c/<l?o ~ 0.1. The red dashed lines are the 
LLs of an infinite graphene sheet which are shifted up (down) 
by Ub = 0.25 eV {-Ub). 



probability density corresponding to the upper states are 
spread out over the dot in both x- and y-directions (see 



Fig. [5jb)) or along the zigzag edges (see Fig. [sjc)). In 
the presence of a p-n junction the electrons confine at 
the p(n)-region when E > Ub and the energy state de- 
crease (increase) with Ub (Fig. [Sjjd)). In contrast with 
Fig. |2j for a p-n junction parallel to the armchair edge 
several states are found for E < Ub- These states, as seen 
in Figs. 5]Je,f), present an interesting behavior: the prob- 
ability densities show a significant localization at the in- 
tersection of the p-n interface and the zigzag edges. That 
can be explained as resulting from the hybridization of 
the zigzag edge states on each side of the p-n junction. 
This remarkable state appears only when the p-n junc- 
tion crosses a zigzag edge. Note that the wavefunction of 
this state consists of both electron and hole components. 



IV. SNAKE STATES: INFLUENCE OF A 
PERPENDICULAR MAGNETIC FIELD 

A. p-n junction 

In the presence of an external magnetic field (see up- 
per panel in Fig. |6] for <i>c/$o = 0.1) the energy spec- 
trum shows anti-crossings for the energies below the gate 
voltage amplitude {E < \Ub\) which is due to the over- 
lap between the quantum Hall (QH) edge states and the 
localized states at the p-n interface (i.e. snake states). 
Because of the smallness of the dot a large magnetic field 
(i.e. B = 800 T for <i>c/$o = 0.1) is required in or- 
der to have a significant influence on the energy levels. 
Nevertheless, as the influence of the magnetic held scales 
with the magnetic flux through the dot area, similar re- 
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FIG. 9: (Color online) (a) Energy levels of a rectangular GQD 
subjected to a p-n junction parallel to the zigzag edges (see 
the inset of Fig. |2| as function of magnetic flux threading 
one carbon hexagon "I>c for Ut = 0.25 eV. The green dashed 
curves are the Landau levels of an infinite graphene sheet 
which are shifted up (down) by Ut {—Ub)- (b) Persistent 
current corresponding to the first electron state (gray solid 
curve, labeled by (I)) and the first hole state (red dashed 
curve, labeled by (2)) as function of external magnetic flux 

suits will be obtained for lower magnetic fields if a larger 
graphene dot is considered. Notice that the number of 
degenerate levels with E = zLUh {E = 0) in the presence 
(absence) of a p-n junction does not change with mag- 
netic field. The red dashed lines in Fig. [6] indicate the 
Landau levels (LLs) of an infinite graphene sheet that are 
shifted up(down) in the presence of an external potential 
Ub{—Ub)- The LLs are given by 

E„ = sgn{n) ^ ± Ut (5) 

where Is = y^h/eB is the magnetic length and n is an 
intcgci^^. Lower panels in Fig. [g] show the probabil- 
ity density corresponding to the states indicated by the 
arrows in the upper panel. Panel (a) shows the confine- 




FIG. 10: (Color online) Probability densities corresponding 
to the points indicated by (a-d) in the enlarged circle in Fig. 
^a). The black line indicates the position of the p-n junction. 




FIG. 11: (Color online) The same as Fig. [ojbut for a rectangu- 
lar GQD subjected to a p-n junction parallel to the armchair 
edges (see the inset of Fig. |4| 
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FIG. 12: (Color online) The position of the oscillations in 
Figs. [9|b) and |ll[b) and the number of flux quanta through 
the surface area S = L xLy 1 LxLy l'2,^ S — - LxLyl\.yi &s func- 
tion of magnetic flux. 
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FIG. 13: (Color online) Probability densities corresponding 
to the points indicated by (a-d) in the enlarged circle of Fig. 
11a). The black line indicates the position of the p-n junc- 



tion. 



ment due to the QH edge states and zigzag edge states 
for Uh — 0. In the presence of a p-n junction and for 
E < Ub, states can arise due to the overlap of the QH 
edge states and the snake states (see panel (b)) or due to 
the overlap with confined states at the p (or n) regions 
(see panel (d)). For E > Ub the carriers form LLs in the 
p (or n) potential regions (see panel (c)). 

The different types of states become more apparent in 
Fig. [7] where we show the current density profile corre- 
sponding to the states shown in the lower panels of Fig. 
[6j The current density vectors are obtained using 

j/^™ = ^[{iplltlmlij^) - (V'ml^LIV'/)] (6) 

where ji^m is the current flowing out of site I into site 
TO. For clarity we show only the current corresponding to 



the A sublattice. Figures [7][b,d) clearly demonstrate the 
presence of snake states at the p-n junction where we have 
clockwise and counterclockwise circling currents, respec- 
tively, in the n and p regions. The current profile also 
reflects the direction of the bonds between the carbon 
atoms and therefore the arrows around the p-n junction 
sometimes point away from the interface. The streamline 
plot in the inset of Fig. ^h) shows the current flow of the 
snake states more clearly. The vector plot in Fig. [Tjja) in- 
dicates the cyclotron orbit of a quantum Hall edge state, 
while Fig. [tJ^c) shows the current profile of a LL state 
that is only very weakly influenced by the p-n junction 
and the edge of the quantum dot. 

The energy levels as function of Ub are shown in Fig. |8j 
in the presence of an external magnetic flux ^c/^o — 0.1 
for the dot with p-n junction along the armchair edges 
(see the inset of Fig. 111). The red solid lines are the LLs 
of a graphene sheet (see Eq. (5)). As in Fig. |6|^a) the 
energy spectrum exhibits different regimes of states; i) 
The regime of QH edge states where Ub < \E\ < E^^ and 
E(-^ is the first LL obtained from Eq. (5). ii) The \E\ < 
Ub and \E\ < E^'" regime where there exist snake states. 
in) The regime of LLs form for > Ef^^. Notice that 
for Ub> {Ub < 0) the LLs form at p (n) region, iv) The 
last regime that can be seen in Fig. [s] (and Fig. [6]ja)) is 
due to the overlap of the edge states and LLs that occurs 
when E^-^^ <\E\<Ub. 

Figure |9][a) displays the energy levels of the system il- 
lustrated in the inset of Fig. [5] as function of magnetic 
flux threading one carbon hexagon <i>j, for Ub = 0.25 eV 
and the same size as Fig. [2] Notice that the zeroth Lan- 
dau level in the absence of a gated voltage is now shifted 
up (down) by +Ub{—Ub)- The green dashed curves are 
LLs of an infinite graphene sheet (given by Eq. (5)). 
The magnetic levels in the GQ D, i.e . the so called Fock 
Darwin levels, approach the LL d^^ ^ ^^ l which are shifted by 
±[/fc. Some of the energy levels approach asymptotically 
the E = zLUb levels. Due to the overlap between the 
QH edge states and the snake states at the p-n junction, 
anti-crossings appear in the energy spectrum. An anti- 
crossing point around ^c/^o = 0.16 is enlarged in Fig. 
|9][a). In Fig. ^h) the persistent current correspond- 
ing to the first electron (solid curve) and the first hole 
(dashed curve) states is shown as function of magnetic 
fiux. The persistent current J is calculated by taking 
the derivative of the corresponding energy levels with re- 
spect to the flux as J{^c) = —dE/d^c- Due to the anti- 
crossings in the energy spectrum, the persistent current 
exhibits an oscillatory behavior with respect to the mag- 
netic flux. The current oscillation due to the snake states 
was recently investigated theoretically for a six-terminal 
graphene nano-ribbon with a p-n junctiorP^. Notice that 
in the presence of the p-n junction the electron QH edge 
states that are shifted down with —Ub overlap with the 
hole QH edge states that are shifted up with Ub (they 
have an opposite circling orbit direction than the elec- 
tronic QH edge states). This hybridize the states in the 
region \E\ < Ub and leads to the current oscillations. 



FIG. 14: (Color online) Current density profile corresponding to the points indicated by (a-c) in the enlarged circle of Fig. 
11a). The black horizental line indicates the position of the p-n junction. 




FIG. 15: (Color online) Energy levels of a rectangular GQD 
with a p-n-p junction parallel to the zigzag edges as function 
of external magnetic flux for the same parameters as Fig. [9] 
The lower inset illustrates the system. The number of p-type 
atoms is twice larger than the number of n-type atoms. The 
green dashed curves are the LLs of an infinite graphene sheet 
which are shifted up (down) by Ut {—Ub)- 



Thus, as the magnetic field is adiabatically increased, at 
each cycle of oscillation the electron becomes predom- 
inantly confined either on the p or the n sides of the 
quantum dot, with the current circulating either clock- 
wise or counterclockwise. 

Figure [TO] shows the electron probability densities cor- 
responding to the points indicated by (a-d) in the en- 
larged region of Fig. |9]^a). Our results indicate that at 
the anti-crossing (panel (b)) the carriers are confined by 
the zigzag edge atoms and the p-n interface which char- 
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FIG. 16: (Color online) (a) Electronic probability density and 
(b) the corresponding current profile for the points indicated 
by (a,b) in Fig. 



15 The black vertical lines indicate the 



position of the p-n junctions. 



acterizes the overlap between the edge and snake states. 
The points corresponding to the energy states that in- 
crease with respect to the magnetic flux around the anti- 
crossing (a,d) are due to states that are confined at the 
p-n junction and the zigzag edges in the n-type region 
(i.e. left side) while panel (c) displays an electron den- 
sity that is confined at the right side of the p-n interface. 

The energy spectrum of the system depicted in the 
inset of Fig. |4]is shown in Fig. [Tl|a) as function of mag- 
netic flux for Ub — 0.25 eV. Since the number of p-type 
and n-type atoms are unequal here (where their mini- 
mum difference is No) the energy levels are not symmetric 
around E = Q. Now the conflnement due to the p-n June- 
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FIG. 17: (Color online) The same as Fig. 15 but for a rectan- 
gular GQD with a p-n-p junction along the armchair edges. 
The system is illustrated in the lower inset. 



tion is along the x-direction which is perpendicular to the 
edge states (caused by the zigzag edges). Therefore the 
energy spectrum exhibits a distinct behavior from that of 
the p-n junction along the zigzag edges. The persistent 
current J corresponding to the energy levels indicated by 
(1) and (2) is shown in Fig. [TI|b) as function of mag- 
netic flux. Notice that the oscillatory behavior is different 
from the results in Fig. [9]^b), now the current amplitude 
decreases smoothly with increasing magnetic flux. The 
position of the oscillations are plotted in Fig 



compared with the flux through the quantum dot (ma- 
genta stars) and half of the quantum dot (blue triangles) . 
Notice that the numerical results are between these two 
curves. It implies that the effective surface area encircled 
by the current is larger than the size of the n or p region. 
The best fit is obtained for flux through a surface area 
S — LrcLy/1.37 (see gray squares). 

Figure [13] shows the electron probability densities cor- 
responding to the indicated points by (a,b,c,d) in the 
enlarged circle (an anti-crossing point around $c/$o = 
0.16) of Fig. Illja). We have a superposition of three 
types of statesTlj zigzag edge states (modified by the p- 
n junction), ii) QH edge states where we have skipping 
orbits, and in) snake states. As seen in the figure the 
overlap of the confined electron in the snake state and 
the edge states leads to a large density at the intersec- 
tion of the p-n junction and the zigzag edges. In contrast 
with the results in Fig. [lO] only half of the zigzag edge 
atoms are contributing to the confinement due to the 
edge states. Therefore the carriers are weakly affected 
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FIG. 18: (Color online) Probability densities corresponding 
to the points indicated by (a-d) in the enlarged region of Fig. 
|17| The black horizontal line indicates the position of the p-n 
junctions. 




,. ,. 



FIG. 19: (Color online) Current density profile corresponding 
to the states shown in Fig. 18 a) and Fig. 18 b). The black 



121 and horizontal lines indicate the position of the p-n junctions. 



by the edge states in comparison with the confinement 
due to the snake states. At the anti-crossing (panel (b)) 
the electrons are mostly confined at the p-n interface and 
along both lengths of the zigzag edges. The correspond- 
ing current profiles of Figs 



13 'a-c) are shown in Fig. 
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Figure[T4{^b) displays the snake states at the p-n interface 
and Figs. [l4]^a,c) show the cyclotron orbit of QH edge 



states respectively at the p and n regions. 



B. p-n-p junction 

Next we investigate the effect of multiple p-n junctions 
where we limit ourselves to the example of two junctions. 
We want to know if there can be any interplay between 
the two junctions, i.e. can states be confined over the 
two junctions? Will there be circling currents between 
the two junctions? 



10 




FIG. 20: (Color online) Energy levels of a rectangular GQD 
with triangular shaped p-n junction as function of external 
magnetic flux for the same parameters as Fig. |9] The lower in- 
set shows schematically the system. The green dashed curves 
are the LLs of an infinite graphene sheet which are shifted up 
(down) by Uh {—Ub). 




FIG. 22: (Color online) Energy levels of a rectangular GQD 
with triangular shaped p-n-p junction as function of exter- 
nal magnetic flux for the same parameters as Fig. [9] The 
lower inset shows the system schematically where the yellow 
regions with zigzag edges indicate the n-type atoms. The 
green dashed curves are the LLs which are shifted up (down) 
by Ub i-Ub). 




FIG. 21: (Color online) Probability densities corresponding 
to the points indicated by (a-d) in the enlarged region of Fig. 
|20| The black line indicates the position of the p-n junction. 



A schematic illustration of a p-n-p junction in a rect- 
angular GQD is depicted in the inset of Fig. [15] where 
the p-n-p junction is parallel to the zigzag edges (along 
y-direction) . The corresponding spectrum in Fig. [Is] ex- 
hibits quite distinct anti-crossings from the case of the 
p-n junction. On the other hand in low magnetic fields 
the hole edge states (i.e. those hole states that decrease 
with respect to the magnetic flux) do not approach the 
zeroth Landau Level (E = —Ub) which is a consequence 
of the fact that the n-type region does not have a bound- 
ary with the zigzag edges. Notice that the hole edge 
states approach E = —Ub in high magnetic fields. The 
electron probability densities for the points indicated by 
(a) and (b) are shown in Figs. [I6|a,b) where the densities 
are spread out mostly along the p-n and n-p interfaces. 
The corresponding current profiles are plotted in Figs. 
|16| ^c,d). Our results show opposite circling currents be- 
tween the two junctions. 

Figure [17] displays the energy spectrum for the p-n- 
p junction parallel to the armchair edges in a rectan- 
gular GQD. The system is depicted in the lower inset. 
Here the n-type region is connected to the zigzag lengths 
which leads to the convergence of the hole edge states 
to i? = —Ub even in low magnetic fields. As in pre- 
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(d) 





FIG. 23: (Color online) Probability densities corresponding to the points indicated by (a-h) in the enlarged region of Fig. 22 
The black lines indicate the position of p-n junctions. 



vious cases for the region \E\ < Ub anti-crossings ap- 
pear in the spectrum. The electron probability densities 
corresponding to the points around the enlarged anti- 
crossing (yellow circle) are shown in Fig. 18 Note that 



the anti-crossing behaviour is qualitative distinct from 
the previous cases shown in Figs. |9|^a) and 11 ^a). Our 
results indicate that at the anti-crossing (a,cj the elec- 
trons are confined due to the edge states corresponding 
to the zigzag atoms in the n-type region and the snake 
states. The points (b,d), correspond to states that ate 
confined at the zigzag edges in the p-type regions and 
near the p-n-p junction. In all points, the probability 
density shows strong peaks at the intersection of the p- 
n junction with the zigzag edges. The density profiles 
have a dumb-bell shape. As before, these localized states 
are associated with the hybridization of the zigzag edge 
states on each side of the junction. However, for non- 
zero magnetic field these states overlap with snake states 
that propagate along the potential interfaces. The cur- 
rent profile corresponding to Figs. [T8|^a,b) is shown in 
Fig. [l9{^a,b) where the counter-circling cyclotron orbits 
in the n and p regions demonstrate the existence of snake 
states at the p-n interfaces. 

V. TRIANGULAR SHAPED P-N JUNCTION 

Next we consider the effect of the gate shape on the 
energy spectrum. Figure [20| displays the energy levels for 
the system illustrated in the lower inset where a triangle- 
shaped gate voltage is assumed for the n-type region. No- 
tice that here we choose an arbitrary direction for the p-n 
junction and it not necessarily matched with the zigzag 
or armchair direction. Since the number of n-type and 



p-type atoms are unequal the electron and hole energy 
levels are not symmetric, i.e. |£'e(*I*c)| \Eh{^c)\- Due 
to the confinement by both edge states (at zigzag edges) 
and snake states (at p-n interface) anti-crossings appear 
in the energy spectrum. An enlargement around one of 
the anti-crossings at <i>c/$o = 0.15 is shown in the in- 
set of Fig. [TT|^b) . The electron probability densities for 
the points around this anti-crossing (a,b,c,d) are shown in 
Fig. 21 Panel (b) shows the electron density at the anti- 



crossing where the electrons confined along the zigzag 
edge and the p-n interface. For the points whose energy 
increases with flux (a,d) the electron is localized along the 
zigzag edge in the p region and snake states are present 
along the p-n junction while the probability density for 
the point (c) is mostly along the zigzag edges in the n 
region and the p-n interface. 

As a last example we investigate the energy spectrum 
of a system consisting of a point contact in a rectan- 
gular GQD (see the lower inset of Fig. 22 1. Recently, 
transport measurements of such a system were carried 
out, and it was found that, due to a chaotic mixing of 
edge channels an unexpected half-integer plateau was ob- 
served in the QH resistivitjP^. The spectrum exhibits 
double anti-crossings between the energy levels in the re- 
gion \E\ < Ub- The enlarged rectangle in Fig. 22 shows 
one of these double-anticrossings around $c/'&o ~ 0.15. 
Figure [23] shows the electron probability densities cor- 
responding to the points indicated by a-h in Fig. [22] 
The energy levels between —Ub and Ub that increase with 
magnetic flux are due to the overlap of the QH edge states 
in the p region and the snake states at the p-n interface 
(Figs. [23][a,b)). Those energy levels that decrease with 
magnetic flux correspond to the zigzag edge states which 
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hybridize with the QH edge states in the n-region and 
to snake states (Figs . [23| [e,g)). Notice that at the anti- 
crossings, see Figs. |23[c,d,f,h), we have an overlap of 
three types of locaHzed states (i.e. QH edge states, snake 
states and zigzag edge states). 

VI. CONCLUDING REMARKS 

We presented numerical results for the energy spec- 
trum and magnetic field dependence of the eigenstates 
of graphene-based quantum dots, on which p-n junctions 
create electron and hole-doped regions. The presence of 
the magnetic field, together with the coupling between 
electron and hole states across the potential barrier due 
to Klein tunneling leads to the appearance of localized 
states at the potential interface, known as snake states. 
These states, which have previously been investigated for 
pn junctions on infinite graphene sheets, can infiuence 
the transport properties of graphene-based nanodevices. 
We have obtained results that show that for the case of 
quantum dots the low energy dynamics of the system 
is dominated by hybridized states that arise due to the 
overlap between quantum Hall edge states and the snake 
states at the p-n junction, with the snake states allowing 
the superposition of quantum Hall edge states at the p 



and n sides of the dot. These states are characterized by 
an energy spectrum that displays an oscillating behavior 
as function of the electrostatic potential and magnetic 
field at the vicinity of the Fermi energy. Furthermore, 
the energy spectrum was shown to depend on the spe- 
cific alignment of the potential interfaces with regard to 
the graphene lattice, as well as on the geometry of the 
gates. The dots were assumed to be defect-free and to 
have perfect zigzag or armchair edges. Future work shall 
concentrate on the effect of edge disorder, impurities and 
defects on the electronic properties of these structures. 
Another aspect that shall be considered is the influence 
of the particular choice of the potential profile and shape 
of the graphene flake on the confined states. 
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